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I.  INTRODUCTION 

An  important  part  of  the  analysis  of  acoustic  waves  in  the  ocean  is 
the  prediction  of  wave  interactions  with  the  bottom.  This  is  a  most 
complex  problem,  not  only  because  of  the  irregular  topography  and 
inhomogeneity  of  the  ocean  bottom,  but  because  the  material  comprising 
the  bottom  is  a  fluid  saturated  porous  material  having  very  complicated 
mechanical  characteristics. 

In  this  report,  a  summary  is  given  of  the  theoretical  models  which 
have  been  developed  for  the  analysis  of  wave  propagation  in  saturated 
sediments.  In  Section  II,  theories  for  wave  propagation  in  the 
unbounded  medium  are  considered.  The  single  component  approach  is 
treated  briefly,  and  then  a  more  extensive  discussion  is  given  of  the 
mixture  theory  due  to  Biot.  In  Section  III,  some  recent  progress  on  the 
reflection  and  refraction  of  waves  at  a  plane  interface  between  water 
and  a  water  saturated  sediment  is  described. 
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II.  WAVES  IN  UNBOUNDED  MEDIA 
A.  The  Equivalent  Continuum 

The  simplest  approach  to  analyzing  waves  in  saturated  sediments  is 
to  model  the  saturated  sediment  as  a  single  continuous  medium.  The 
mechanical  properties  of  the  continuous  medium  are  chosen  so  that  its 
behavior  is  equivalent,  in  a  necessarily  limited  sense,  to  that  of  the 
saturated  sediment.  Since  waves  in  saturated  sediments  in  general 
exhibit  dispersion  and  attenuation,  the  equivalent  continuous  medium 
must  be  viscoelastic. 

1-5 

Models  of  this  type  have  been  discussed  by  Hamilton.  He  suggests 

the  use  of  the  particularly  simple  linearly  viscoelastic  model  in  which 

the  Lame  constants  X  and  /a  of  an  isotropic,  linearly  elastic  material 

are  replaced  by  complex  moduli  X+ix'  ,  •  The  moduli  X,  X',  , 

and  fji''  are  assumed  to  be  frequency  dependent  and  are  to  be  determined 

from  experiment.  This  type  of  model  is  equivalent  to  a  Voigt  material.® 
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Hamilton  cites  developments  of  the  model  by  Ferry  and  White. 

Such  a  purely  empirical  approach  has  substantial  advantages  and  dis¬ 
advantages.  The  primary  advantage  is  its  simplicity.  Once  the 
coefficients  in  the  theory  have  been  evaluated  as  functions  of 
frequency,  then  the  large  body  of  literature  and  many  existing  computer 
programs  on  wave  propagation  in  viscoelastic  materials  can  be  used  to 
obtain  solutions  to  specific  problems.  On  the  other  hand,  since  the 
model  does  not  explicitly  incorporate  the  physical  properties  of  the 
saturated  sediment,  the  coefficients  in  the  theory  must  be  measured  anew 
for  each  change  in  the  properties  of  the  material.  More  important,  the 
theory  does  not  have  the  ability  to  predict  the  effects  on  wave 
propagation  of  changes  in  the  material  properties. 
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B.  The  Biot  Theor 


In  order  to  introduce  the  physical  properties  of  the  saturated 
sediment,  the  material  can  be  modeled  as  a  binary  mixture  of  fluid  and 
solid  constituents.  By  retaining  the  individual  identities  of  the  fluid 
and  solid,  it  is  possible  to  express  at  least  part  of  the  coefficients 
in  terms  of  the  physical  properties  of  the  two  constituents. 


A  theory  of  this  type  was  developed  for  application  to  fluid 

Q  IP 

saturated,  porous  elastic  media  by  Biot.  While  he  derived  his 

equations  on  a  somewhat  intuitive  basis,  many  subsequent  studies  have 
confirmed  that  they  are  consistent  with  the  principles  of  continuum 
mechanics  (see,  for  example,  Bowen^^  and  Bedford  and  Drumheller^^). 


Biot's  equations  can  be  written 


(l-(i.)  Pg  U(s)n,  -  -"’(‘J(s)m‘^f)m^ 


^  ®(s)kk,m  ®(s)mk,k  ^  ^  ®(f)kk,m  '  ^^‘^(s)m"“(f)m 


■  ’  (1) 


•^Pf  ^'f)m  "  "’(“(s)m'‘^(f)m) 

"  Q  ®(s)kk,m  "  "(f)kk,m  ^  ^(^s)m-^f)m)  ’  (2) 

where  is  the  porosity,  and  p^  are  the  mass  densities  of  the 
solid  and  fluid,  and  are  the  displacement  vectors  of  the 

solid  and  fluid,  the  notation  ,m  denotes  partial  differentiation  with 
respect  to  the  coordinate  x^^^,  and  are  the  linear  strains 

of  the  solid  and  fluid  [e(s)mk*l/2(u(s)^^k^U(s)k,m)]  •  A,  Q. 
c,  and  b  are  constitutive  coefficients. 
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Equations  (1)  and  (2)  are  the  equations  of  motion  for  the  solid  and 
fluid  constituents,  respectively.  The  first  term  on  the  right  of  each 
equation  is  a  virtual  mass  term  which  is  linear  in  the  relative 
acceleration  of  the  constituents.  The  last  term  on  the  right  is  a  drag 
term  which  is  linear  in  the  relative  velocity  of  the  constituents.  By 
using  the  solution  for  an  oscillating  cylinder  containing  a  viscous 
fluid,  Biot  was  able  to  evaluate  the  coefficient  of  the  drag  term  b  as  a 
function  of  frequency.^® 

In  a  further  important  development,  Biot  and  Willis^®  showed  that 
the  coefficients  A,  Q,  and  R  could  be  expressed  in  terms  of  the  bulk 
modulus  of  the  solid  material,  the  bulk  modulus  of  the  fluid,  and  the 
bulk  and  shear  moduli  of  the  drained  porous  solid  and 

The  Biot  theory  was  first  applied  to  marine  sediments  by  Stoll  and 
Bryan^®  and  Stoll. In  order  to  account  for  dissipative  effects 
associated  with  motion  of  the  granular  matrix  material,  these  authors 
assumed  that  the  moduli  K|j  and  were  viscoelastic.  Based  upon 
experimental  data  on  dry  granular  media,  it  was  assumed  that  K.  and  m. 
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were  complex  constants.  It  was  shown  by  Stoll  that  the  theory  could 
correctly  predict  the  variation  with  frequency  of  the  attenuation  of 
compressional  waves  in  saturated  sands. 

The  work  of  Stoll  and  Bryan  was  extended  by  Hovem  and  Ingram,^®  who 
showed  that  Biot's  technique  for  evaluating  the  drag  coefficient  b  could 
be  used  to  evaluate  the  virtual  mass  coefficient  c  as  well.  They  were 
also  able  to  obtain  a  more  explicit  expression  for  b  in  the  case  of  a 
sediment  consisting  of  particles  of  uniform  size.  They  then  presented 
very  favorable  comparisons  of  the  theory  with  measurements  of 
compressional  wave  attenuation  in  saturated  sands.  In  addition,  they 
presented  new  experimental  data  on  compressional  waves  in  a  model 
sediment  consisting  of  spherical  glass  beads  saturated  by  water,  and 
showed  that  the  theory  correctly  predicted  both  the  attenuation  and 
phase  velocity  of  compressional  waves. 


Thus  the  extended  Biot  theory  has  been  shown  to  give  quite  accurate 
predictions  in  comparison  with  direct  measurements  of  both  phase  j 

velocity  and  attenuation  in  saturated  sediments. 

Stoll^^  has  used  the  Biot  theory  to  show  that  the  dominant  loss 
mechanism  at  high  frequencies  is  the  drag  between  the  solid  and  fluid 
constituents,  while  the  dominant  mechanism  at  low  frequencies  is 
dissipation  in  the  viscoelastic  granular  matrix. 

Because  the  scale  of  measurements  of  phase  velocity  and  attenuation 

is  necessarily  small  in  a  laboratory  setting,  the  data  are  in  the 

kilohertz  frequency  range  and  above.  Such  measurements  cannot  be  used 

PI 

to  verify  the  theory  at  lower  frequencies.  Stoln  has  obtained  some 
low  frequency  data  on  the  attenuation  of  shear  waves  in  saturated 
sediments  by  using  a  resonant  column  technique,  but  he  gave  only  * 

qualitative  comparisons  with  the  theory. 

There  is  a  pressing  need  for  additional  data  at  low  frequencies  and 
for  data  which  can  shed  light  on  the  viscoelastic  properties  of  the 
granular  matrix. 

22 

Recently,  Shirley  et  al.  presented  direct  measurements  of  phase 
velocity  and  attenuation  in  a  glass  bead  sediment  saturated  by  a  mixture 
of  water  and  glycerine.  By  varying  the  proportion  of  glycerine,  the 
viscosity  of  the  liquid  could  be  changed.  In  comparison  with  the 
extended  Biot  theory,  the  data  predicted  a  substantially  greater 
increase  in  attenuation  with  increasing  viscosity.  The  measurements 
were  made  at  a  single  frequency.  Measurements  made  for  a  range  of 
frequencies  would  greatly  aid  in  characterizing  the  viscoelastic 
response  of  the  granular  matrix.  This  in  turn  could  conceivably  permit 
extrapolation  of  the  theory  to  low  frequency. 
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.  Reflection  at  an  Interface 


The  ultimate  objective  of  efforts  to  model  wave  propagation  in 
saturated  sediments  is  the  analysis  of  wave  interactions  with  the  ocean 
bottom.  Although  much  work  remains  to  be  done  on  propagation  in  the 
unbounded  medium,  important  progress  on  the  interaction  problem  has 
recently  been  made  by  Stoll  and  Kan.^^  They  considered  a  liquid  half 
space  above  a  half  space  of  saturated  sediment,  and  solved  the  problem 
of  the  reflection  and  refraction  of  plane  waves  at  the  interface  using 
the  extended  Biot  theory.  The  analysis  was  based  on  the  earlier  work  on 

OA 

saturated  porous  media  by  Deresiewicz  and  Rice. 

In  contrast  to  the  results  for  elastic  media,  they  showed  that  the 
reflection  coefficient  exhibited  strong  frequency  dependence,  and  was 
also  a  function  of  the  assumed  value  of  permeability  of  the  sediment. 
The  results  obviously  have  important  implications  for  the  modeling  of 
bottom  interactions. 

The  analysis  of  Stoll  and  Kan  will  have  to  be  extended  before 

realistic  estimates  of  bottom  loss  can  be  achieved.  They  assumed  that 

the  sediment  properties  were  homogeneous  and  isotropic.  It  is,  of 

course,  well  known  that  sediment  properties  are  strongly  depth 

5  19 

dependent.  ’  There  is  also  some  evidence  that  sediments  are 
25  26 

anisotropic.  ’  Both  of  these  factors  would  influence  the  reflection 
and  refraction  of  waves. 
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